Detection of TeV Gamma Rays from the BL Lac Object 
1ES1959+650 with the Whipple 10m Telescope 

J. Holder 1 , I.H. Bond 1 , P.J. Boyle 2 , S.M. Bradbury 1 , J.H. Buckley 3 , D.A. Carter-Lewis 4 , 

W. Cui 5 , C. Dowdall 6 , C. Duke 7 , I. de la Calle Perez 1 , A. Falcone 5 , D.J. Fegan 8 , 
S.J. Fegan 9 ' 10 , J.P Finley 5 , L. Fortson 2 , J.A. Gaidos 5 , K. Gibbs 9 , S. Gammell 8 , J. Hall 11 , 
T.A. Hall 12 , A.M. Hillas 1 , D. Horan 9 , M. Jordan 3 , M. Kertzman 13 , D. Kieda 11 , J. Kildea 8 , 
J. Knapp 1 , K. Kosack 3 , H. Krawczynski 3 , F. Krennrich 4 , S. LeBohec 4 , E.T. Linton 2 , 
J. Lloyd-Evans 1 , P. Moriarty 6 , D. Mliller 2 , T.N. Nagai 11 , R. Ong 14 , M. Page 7 , 
R. Pallassini 1 , D. Petry 4 ' 16 , B. Power- Mooney 8 , J. Quinn 8 , P. Rebillot 3 , PT. Reynolds 15 , 
H.J. Rose 1 , M. Schroedter 9 - 10 , G.H. Sembroski 5 , S.P Swordy 2 , V.V. Vassiliev 11 , 
S.P Wakely 2 , G. Walker 11 , T.C. Weekes 9 



ABSTRACT 



We present the first strong detection of very high energy 7-rays from the 
close [z = 0.048) X-ray selected BL Lacertae object 1ES1959+650. Observations 
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were made with the Whipple 10m telescope on Mt. Hopkins, Arizona, using the 
atmospheric Cherenkov imaging technique. The flux between May and July 2002 
was highly variable, with a mean of 0.64 ± 0.03 times the steady flux from the 
Crab Nebula and reaching a maximum of five Crab, with variability on timescales 
as short as seven hours. 

Subject headings: BL Lacertae objects: individual (1ES1959+650) — gamma 
rays: observations 

1. Introduction 

Blazars, consisting of flat spectrum radio quasars (FSRQs) and BL Lacertae objects 
(BL Lacs), are a relatively rare type of active galaxy characterized by extremely variable, 
non-thermal spectral energy distributions (SEDs). The emission extends from radio to 7- 
ray frequencies and is believed to be produced by a highly relativistic plasma jet closely 
orientated with the line of sight to the galaxy (Urry and Padovani 1995; Blandford and Rees 
1978). In a vF v representation, the SEDs display two peaks. The lower energy peak is usually 
attributed to synchrotron radiation from relativistic electrons, while the higher energy peak 
is thought to be due to inverse Compton scattering by these same electrons. The photons 
involved in the scattering process may be either synchrotron photons produced in the jet, in 
the synchrotron self-Compton (SSC) model, or some external population, such as photons 
emitted from an accretion disc (Dermer, Schlickeiser and Mastichiadis 1992) or reflected 
from emission line clouds (Sikora, Begelman and Rees 1994). Alternative explanations for 
the high energy component invoke pair cascades triggered via pion and pair photoproduction 
from high energy protons in the jet (Mannheim 1993, 1998), or proton synchrotron radiation 
(Aharonian 2000; Miicke and Protheroe 2001). 

The EGRET detector on board the CGRO detected in excess of 66 blazars (Hartman 
et al. 1999) above 100 MeV; however, only a few have been detected at higher energies 
(> 300 GeV) by ground based atmospheric Cherenkov telescopes. Markarian 421 (Punch 
et al. 1992) and Markarian 501 (Quinn et al. 1996) were the first extragalactic TeV sources 
to be detected and the behaviour of these two nearby (z = 0.031 and z = 0.034 respectively) 
objects has been closely monitored and studied. Flux variability over two orders of magnitude 
has been observed, with doubling timescales as short as ~ 15 minutes (Gaidos et al. 1996), 
placing strong constraints on the size of the emission region. The 7-ray flux is correlated 
with the X-ray flux (Buckley et al. 1996; Catanese et al. 1997a; Krawczynski et al. 1997) 
and, in the case of Markarian 421, with the spectral power law index in the 7-ray region 
(Krennrich et al. 2002; Aharonian et al. 2002a). Unconfirmed detections of 1ES 2344+514 
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(Catanese et al. 1998) and PKS 2155-304 (Chadwick et al. 1999) have also been reported 
and, most recently, the detection of very weak emission from H1426+428 has been confirmed 
by three groups (Horan et al. 2002; Aharonian et al. 2002b; Djannati-Atai et al. 2002). This 
last object is of particular interest since it is relatively distant (z = 0.129) and hence the 
observed flux is expected to be strongly attenuated by absorption via pair production on 
the extragalactic infra-red background light (EBL) (Gould and Schreder 1967; Stecker, De 
Jager and Salamon 1992). 

All of the TeV blazars can be classified as high frequency peak BL Lacs (HBLs) (Padovani 
and Giommi 1995) on the basis of the location of the lower energy peak in their SEDs, and 
predictions of TeV emission from other nearby sources of this type have driven observations 
by atmospheric Cherenkov detectors. 1ES1959+650 (z = 0.048) was first suggested as a good 
TeV candidate of this type by Stecker, De Jager and Salamon (1996), who used simple scaling 
arguments to predict it as the third strongest TeV blazar, after Markarian 421 and Markar- 
ian 501, with a flux prediction of 1.9 x 10~ n cm~ 2 s _1 above 300 GeV. More recently Costa- 
mante and Ghisellini (2002) have predicted fluxes above 300 GeV of 7.5 x 10~ n cm~ 2 s _1 , 
based on a simple phenomenological parameterization of the SED adapted from Fossati et al. 
(1998), and 0.03 x 10~ n cm~ 2 s _1 using a homogeneous, one- zone synchrotron self-Compton 
model (Ghisellini, Celotti and Costamante 2002). 

In fact, 1ES1959+650 is not a particularly extreme example of an HBL. Measurements 
with BeppoSAX in 1997 (Beckmann et al. 2002) placed the lower energy SED peak at 
10 15 Hz (4 eV), as compared to typically ~ 1 keV for Markarian 421 (Maraschi et al. 1999) 
and ~ 100 keV for Markarian 501 during its 1997 flare state (Pian et al. 1998), although 
for both of these sources the position of the peak is known to vary widely depending upon 
the flux level. Giebels et al. (2002) report on X-ray observations of 1ES1959+650 obtained 
with the USA(1-16 keV) and RXTE (2-16 keV) missions during 2000, showing threefold 
increases in the X-ray flux on a timescale of a few days. The flux increase is correlated 
with spectral hardening, indicative of a shift of the lower energy peak of the SED towards 
higher frequencies. The source is also unusual among HBLs for its strong rapid optical 
variability. Observations by Villata et al. (2000) showed rapid flickering, including a decrease 
of 0.28 magnitudes in four days. 

Prior to May 2002, only tentative evidence had been presented for TeV emission from 
1ES1959+650. A detection with a statistical significance of 3.9 a was reported by Nishiyama 
et al. (1999) based on 57 hours of observations with the Utah Seven Telescope Array. More 
recently, Konopelko et al. (2002) reported a preliminary detection at ~ 5 a for the HEGRA 
Cherenkov telescope system. Previous observations of 1ES1959+650 with the Whipple 10m 
telescope produced an upper limit at a flux level of 1.3 x 10~ n cm -2 s -1 above 350 GeV 
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(Catanese et al. 1997b). 

Observations of 1ES1959+650 during May- July 2002 with the Whipple 10m telescope 
resulted in a clear detection of TeV 7-ray emission from this object (Dowdall, Moriarty and 
Kosack 2002) which was rapidly confirmed by the HEGRA experiment (Horns and Konopelko 
2002). We report here the results of the Whipple 10m observations. 

2. Observations and Data Analysis 

The configuration of the Whipple 10m 7-ray telescope during these observations is 
described in detail in Finley et al. (2001). Briefly, the telescope consists of a 10 m reflector 
and a 490 pixel photomultiplier tube (PMT) camera. For the purposes of this analysis 
only the high resolution (0.12° spacing) central 379 PMT pixels have been used. Cherenkov 
images are recorded and parameterized according to Hillas (1985) and then 7-ray like images 
are selected using the "supercuts" criteria (Reynolds et al. 1993) optimized on recent data 
from the Crab Nebula, the standard candle of TeV 7-ray astronomy (Weekes et al. 1989). 
Following a realignment of the optical system in February 2002, observations of the Crab 
Nebula showed that the telescope was at its most sensitive since the installation of the 
current camera, with one hour of Crab observations producing a 6a detection; however, 
between February and July 2002 a decrease in the relative efficiency of the telescope of 
~ 30% was measured by examining the response to the cosmic ray background. The cause 
of the effect is still under investigation - one explanation may be that it is due to increased 
atmospheric absorption caused by large forest fires in the region during this period. 

Observations of 1ES1959+650 were made during moonless periods between May 16 th 
and July 8 th 2002. The total dataset consists of 39.3 hours of on-source data, together with 
7.6 hours of off-source data for background comparison. For observations from Mt. Hopkins 
(latitude 31°57.6'N) 1ES1959+650 culminates at a zenith angle (9 Z ) of 33.5° and so the data 
were necessarily taken at large 9 Z , between 53.5° and 33.5°. After accounting for the zenith 
angle and reduced telescope efficiency, we estimate the energy threshold {E t h res h) for the 
majority of the observations to be ~ 600 GeV . E t hresh is defined here as the energy of 
the peak differential 7-ray flux for a source with the same spectrum as the Crab Nebula. 
Figure 1 shows the distribution of the Cherenkov image orientation angle alpha for the full 
on-source dataset, together with the distribution for the off-source dataset scaled such that 
the number of events in the region a > 30° is the same. The average rate of the excess in 
the on-source 7-ray region at a < 15° is 1.08 ± 0.05 7 min -1 , corresponding to a detection 
of greater than 20 a. 
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3. Flux Variability 

Figure 2 shows the daily averaged rates for 1ES1959+650 for all of the observations. 
The rates are expressed in multiples of the steady Crab Nebula flux and have been corrected 
in order to account for the varying zenith angle and changes in telescope efficiency according 
to the procedure of LeBohec and Holder (2002). This has been tested using Crab Nebula 
data taken over a wide range of zenith angles and atmospheric conditions. The correction 
factor varies from run to run and has a mean value for this dataset of 2.0 ± 20%. Strong 
night-to-night variability is evident; the largest change in rate, between MJD 52428 and 
MJD 52429, corresponds to a doubling time of 7 hours - shorter than has ever been observed 
in other wavebands for this source. The mean flux over all observations was 0.64 ±0.03 Crab 
units. 

Figure 3 shows the rate in 5 minute bins for two nights, May 17 th and June 4 th , during 
which the source was most active. The statistical evidence for variability within these nights 
is given by the \ 2 probabilities of constant emission = 1% and 8% respectively. We conclude 
that there is no strong evidence for flux variability on this timescale. 

4. Discussion 

The detection of TeV 7-ray emission from 1ES 1959+650 adds another member to the 
class of TeV blazars, all of which are close BL Lac objects having a low bolometric luminosity 
and the peaks in their SEDs at high frequency. The rapid flux variability, and the fact that 
1 ESI 959+650 has not been detected during previous observations, indicates that the source 
was in an unusual flaring state during these observations. The flux level was at times orders of 
magnitude above the most recent model predictions. Throughout the period of the Whipple 
observations measurements in the 2 — 12 keV region by the All-Sky Monitor (ASM) on board 
RXTE have shown 1ES1959+650 to be active and variable, with daily average fluxes reaching 
~ 20 mCrab in May and July. Target of opportunity observations with the RXTE small 
field of view instruments were triggered following the 7-ray detection and will be reported 
on elsewhere (Krawczynski, private communication). The rapid flux variability observed at 
TeV energies implies a small emission region in the jet with a high Doppler factor (Mattox 
et al. 1993; Madejski 1996; Buckley et al. 1996). The contemporaneous X-ray and TeV 7-ray 
data will allow us to constrain the jet parameters when modelling the emission processes. 

Analysis of the Whipple data is ongoing, but attempts to reconstruct the source spec- 
trum have been hampered by the effects of the decreased telescope efficiency during the 
observation period. The HEGRA collaboration measure a rather steep spectrum (spectral 
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index a = 3.2 ±0.3) for observations prior to 2002, while the spectrum during the flaring pe- 
riod exhibits pronounced curvature and deviates significantly from the spectrum seen during 
the quiescent state (Horns et al. 2002). The majority of models for the EBL lead to predic- 
tions of a cut-off in the 7-ray region beginning below ~ 10 TeV for a source at z — 0.048 
(Primack (2001); but see Vassiliev (2000) for a discussion of an EBL model which does not 
produce a distinct feature in the observed spectrum). Deviations from a pure power law 
have now been resolved in the spectra of both Markarian 421 (Krennrich et al. 2001; Piron 
et al. 2001; Aharonian et al. 2002a) and Markarian 501 (Samuelson et al. 1998; Aharonian 
et al. 1999; Djannati-Atai et al. 1999) and the spectrum of the most distant TeV blazar, 
H1426+428, is measured to be very steep (a = 3.50 ± 0.35(stat) ± 0.05(syst)) (Petry et al. 
2002), but it is not yet clear whether these features are due to absorption on the EBL or are 
intrinsic to the sources. Clearly, further observations and spectral analysis of 1ES1959+650 
may help to resolve this question. 

The VERITAS Collaboration is supported by the U.S. Dept. of Energy, N.S.F., the 
Smithsonian Institution, PPARC (U.K.) and Enterprise Ireland. 
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Fig. 1. — Distribution of alpha for all on and off-source events (the off-source distribution 
has been scaled so as to match the number of events in the a > 30° region). Alpha is defined 
as the angle between the major axis of an elliptical 7-ray image and the source location in 
the camera. 



-10- 




52410 52420 52430 52440 52450 



MJD 

Fig. 2. — The daily average 7-ray rates for 1ES1959+650 during 2002. The rates have been 
corrected for zenith angle of observation and relative telescope efficiency as described in the 
text. The dashed line indicates the most rapid rate change, corresponding to a doubling 
time of 7 hours. 



-li- 



es 

•8 
u 

CD 
■<— i 




10 10.5 11 11.5 

UT Time (hours) 




UT Time (hours) 



Fig. 3. — The 7-ray rates for the two nights showing the most activity (5 minute binning). 
The rates have been corrected for zenith angle of observation and relative telescope efficiency 
as described in the text. 



